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I. HISTORY AND INTRODUCTION 

This thesis is the third in a series of investigations 
into a proposal put forth by R. Lauer of Naval Underwater 
Systems Center, New London, CT , in 1979. In his memorandum 
[Ref, 1], lauer describes a "new" way to analyze sound 
propagation in the ocean and two applications of the method, 
which he named the " Wavenumber Technique,' 1 or F.T. In his 
proposal, Lauer stated that a source of sound in the ocean 
could be pinpointed in both depth and range by a single 
omnidirectional hydrophone provided the source is generating 
a continuous wave tone. 

The W. T. was initially described by F. DiNapoli [Ref. 2] 
as an intermediary step in his development of a speedier 
computer algorithm used in analyzing sound propagation as a 
function of range. DiNapoli and Lauer proposed converting 
the pressure as a function of range ^{r) to the pressure 
spectrum as a function of the wavenumber P (k) . The conver- 
sion is accomplished by taking a weighted Fourier transform 
of the pressure field ^(r). Once the pressure spectrum is 
obtained, analysis of the pressure field can be accomplished 
in a manner analogous to that presently used in signal 
processing. 

Lauer proposed two uses for his "Wavenumber Technique:" 

1. determination of the depth of the acoustic source, an 
ability which has obvious tactical applic ations ; 

2. use of the W. 1. to evaluate the accuracy of existing, 

and future, acoustic models such as P.E. and 

F.A.C.T., by breaking the generated pressure field 
down into ray path "families," such as bottom-bounce, 
ref racted-sur face- ref lected, and surf ace- ducted fami- 
lies, thus giving a quantitative read-out of what 



proportion of the total acoustic energy is being 
channelled through the various ray paths. 

For both uses, a knowledge of the acoustic environment is 
necessary. 

It is the intended purpose of this thesis to investigate 
the validity of the first application of Lauer's wavenumber 
technigue by using it in conjunction with an acoustic 
phenomenon that is well-known and for which accurate results 
can be calculated with precision; this phenomenon is the 
ILoyds Hirror for an Acoustic Doublet. Using several 
source/receiver combinations, this investigator intends to 
compare actual source depth with that predicted by the W.T. 
If Lauer’s technigue is cogent, the source depths calculated 
via the two methods should be equal, or very close within a 
statistically acceptable degree. 

As mentioned in the beginning of this thesis, two prior 
inves tigat ors, B. Stamey and J. Blanchard, looked into 
applying the W. T. to determining source depth. Stamey’ s 

investigation [Ref. 3] utilized a parabolic equation 

computer model developed by H. K. Brock, the Split Step Fast 
Fourier Transform, or SSFFT, to generate the pressure 
fields. He qualitatively compared the model to a Normal 
Hoae model and a P. 2. Finite Difference model, both range 
dependent, and to DiNapoli’s Fast Field Program model, a 
range independent model. Stamey concluded that the W.T. 
showed promise as an analysis tool and that further investi- 
gations were warranted. 

J. Blanchard carried Stamey ’s researches one step 
further [Ref. 4]. Using two Parabolic Equation computer 
algorithms, 3rock’s Split Step Fast Fourier Transform and 
Jaeger's Implicit Finite Difference, as his pressure field 
generators, he examined the use of the spacing between nulls 
of the pressure spectrum to determine source depth. His 
results are interesting and support the need for further 
invest igation. 
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satisfied with 



The author was not 



Stanley's and 



Blanchard's findings. 3oth individuals acknowledged short- 
comings in their respective studies, especially with regard 
to the mutual presence of the "U-shaped phenomenon" encoun- 



however, that their theses were preliminary studies only and 
were produced within a highly restricted time frame. Also, 
they used existing computer algorithms specifically designed 
to approximate acoustic propagation in a velocity- variant 
medium. However, the environmental settings used by both 
gentlemen were that of the Lloyds Mirror for an Acoustic 
Doublet which requires a constant sound speed throughout the 
water column; thus, the approximations, assumptions and 
"fudge factors" used by these models make their results for 
an isospeed medium highly artificial and somewhat suspect. 
Fortunately, however, there exists for the Lloyds Mirror a 
simple, geometric solution specifically designed for an 
isospeed environment [Ref. 5]. Since the pressure fields 
for varying source and receiver depths and, therefore, the 
corresponding pressure spectra, can be precisely calculated, 
it was thought this model would be a good check on the oper- 
ational applicability of Lauer's proposed technique. 

Secondary considerations of this study were to examine 
what effect, if any, the introductions of, firstly, surface 
waves and, secondly, noise would have on the pressure spec- 
trum. The first objective was simulated by use of 
mathematical formulae given in reference 5; for the second 
objective, it was thought that complex noise, random in both 
amplitude and phase and similar to a Rician distribution 
[Ref. 6: p. 189], would provide a simple but reasonable 
approximation of an ocean noise field. This investigator 
intends to gradually intensify the surface waves and the 
noise field, separately, until the original pressure spec- 
trum is no longer recognizable. In this manner, it should 



terea in the pressure spectrum 



One should remember. 



16 



be possible to make a qualitative assessment of their 
respective debilitating effects under the carefully 
controlled conditions found in the Lloyds Mirror phenomenon. 

Section II sets forth in more detail the theoretical 
development of each point described in this introduction. 
Section III summarizes the investigator’s results and the 
conclusions drawn from those results. Section IV contains a 
listing of the computer algorithms utilized. 
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II. T3E0RY 

This chapter presents the mathematical basis of the 
wavenumber technique. 

A. THE LLOYDS HIRROB PHEHOMENON 

A diverging mcnof reguency spherical pressure wave 

[Bef. 5: p.112] can be written in complex form as 



,p(R,f) -£(R)e 



-iut = A e i(kR -u f) 

R 



(egn 2.1) 



where ,P (B) is the spatial factor. Referring to Figure 2.1, 
R, is the range from the source to the receiver and B, is 
the range from the image to the receiver. For convenience 
in all tha't follows, we shall set A to unit magnitude. The 
following equations apply to the acoustic waves propagating 
directly out from the source P 5 and from the apparent image 
P ; (which is actually the wave reflected from the surface): 

£/R,t) = £ s (Rle- lul = -fee 



£-JR,t) =£t(R) e -‘ ut - ^e"“ r " 



t) 



(ean 2.2) 



where P , P. , and E» are complex functions of horizontal 
range r^vertical depth Z and time t; the minus sign in the 
equation for P; is derived from the surface reflection coef- 

Inspection of 



ficient, \jjp = -1, 



for a smooth surface. 



Figure 2.1 reveals that 



Rj = V(z / .-z s ) 2 + ^ 2 ^ + Zs ^ + r ? 



1 8 



and A r - 



R-R, 



1.2 



And so the total field can be written as 



- P(R)c 



-iU)t 



J(kR r Ut) 

~~R, 



J(kR r (x)t) 

'~Ro 



(egn 2.3) 



Since P. and P; both have the same time factor, 

<s *5 /vO ’ 

exp ( — Lco*t) ^ can retain just the spatial factors and equa- 

tions 2.3 reduce to 

e ikR, e ikR 2 " 

R, R 2 . 

(eqn 2.4) 



P(R) = 



equation 2.4 is a form of the complex pressure as a function 
of range used in the computer algorithm shown in Appendix I. 

Inspection of Figure 2. 1 reveals that, for R >> Z s and 
Q very small, 

r^R 

or 

sin 9 -jj- 



or 



A '~7T 



(eqn 2.5) 



and the pressure amplitude can be approximated by 



F(R) 





(eqn 2.6) 



Looking at just the formula for the pressure amplitude, 
ve can see that as 
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nn 3 n-0 } l } 2 ? 3 r . . 



(egn 2.7) 




the pressure amplitude goes to zero, producing the classic 
1 P ( H ) 1 vs E curve shown in Figure 2.2 



B. THE BEIATIOHSHIP EETEEEN K, y , AND p 

In his description of the Wavenumber Technique [Ref. 1: 
p. 5-6], Lauer utilizes the horizontal and vertical compo- 
nents of the wavenumber, k. The general relationship amoung 
these three terms is illustrated in Figure 2.3 and can be 
written mathematically as 

k = the wavenumber = 2vf/c 

J - the horizontal component of k = k cos (p 
P = the vertical component of k = k sin (fi 
or 



k •- •)/ y 2 +P 



2 



(egn 2 .8) 



C. THE WAVENUMBER TECHNIQUE AND THE LLOYDS MIRROR 

Given a point source in free space, the monof reguencv 
pressure field at the receiver, ,P(R) » can be expressed as a 
spherical wave (time factor suppressed) , 
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£(R) - £ 



ikR 

R 



where, from Figure 2.1, 



(egn 2.9) 



or 



R 2 = r 2 + Z r 2 , for Z r > 



R 2 = r 2 + Z g 2 , for Z r < Z g 



(egn 2. 10) 



In integral form, P (E) can be written as [Ref. 7: p. 127] 



P(R) = £ 



ikR f 00 J e (7r) e ti P\ Z r~ Z s\ 



R 



n 



dr 



(egn 2.11) 



which is taken from the Fourier-Bessel Transform Pair 
[Ref. 7: p = 126 ], 

fjr) = r c °gf ?; Jjyr) ydy 



Q(r) = J 0 r yr) r dr 



(egn 2. 12) 



where f,(r) represents the acoustic pressure function in the 
range domain, and < a,{ / y) is the acoustic pressure spectrum in 
the wavenumber domain. The sign of the exponential function 
in the integral is based on which of Z 5 and Z r is greater. 
For the case of receiver depth being greater than source 
depth so that waves from both image and source are travel- 
ling downward at the receiver depth, the total pressure at 
the receiver can be expressed as 



E(R) = 



■CO e ip(2 r -Z s ) _ e ip(2 r + Z s ) 

Tp 



JJyr)ydy 



(egn 2. 13) 



Use of Euler's Identity reduces Eguation 2.13 to 
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(egn 2.14) 



£(R) = ~2 f — - - s ' n( P Z s) JJ y r j y cy 

Jq p 

Using the relationship between the Hankel functions [Ref. 5: 
p. 449], 

H ( J } (yr) = JJyr) + iYJyr) 

and 

( 2 ) 

H 9 (yrj = JJyr) - YJyr) 

(egn 2. 15) 

one can re-vrite the Eessel function in equation 2.14 as 

JJyr) = 4; [h/ ( o n (yr) + ri[ 2 \yr)\ (egn 2. 16) 

so that 

2 a. ( y ) = l l<i'>H , ."(yr)rdr * \Jt(r)H',\r)rdr (egn 2 _ 17) 

Letting r’ = -r and looking at the second term on the right 
hand side of equation 2.17, 

l JJ-r J ) h{ 2) (-yr') r'<jr J = J^MH 9 n (yr)r dr (egn 2.18) 
Now, assuming that f (r) = f (-r) , then 

2ajy) = j'JLfr) H['\yr) r dr (egn 2.19) 

Therefore, 
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(eqn 2. 20) 



r co , ,(l), , - 

o (R) = - f jHoJirL jP z r 



X. 



J-CO (2 



sin(pz s )dy 



For yx >2 71 , the asymptotic approximation of th< 

Eankel function for large argument can be used: 



H ( J } (yr) =» 



i(yr^) 



nyr 



(egn 2.21) 



assuming |r (r) | goes to zero faster than (In r) , 



r Z r 00 



e'@ z ' sin(pz s ) 



e‘^ r dy (egn 2.22) 



Let the terms inside the brackets in equation 2.22 be 

defined as g('V) , then ^ 

^ / ' nr r co 



D 

-C 



riJ i 7 r r 

- - V 777 e' T v y 



and 



f g(y)e‘ yr dy = - JZ^E7/7“ £Y/?; (e qn 2. 

7 -CO ^ 
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or, using equation 2.12, 

£fr; = ' 



7T” 

/‘IT 






* OO 



y “F" ^7T / P(R)z ,yr dr 



-co 



(eqn 2.24) 



Note that equation 2.24 has the form of the Fourier trans- 
form cf the pressure function if we define 



Ur) = -fr~ £(R) 



( eqn 2. 25) 
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so that vr p (E ) and ^('y) are Fourier Transform pairs. 
This equation can be easily evaluated [Ref. 8], 



3Jy) 



7H 

/ T" 






\^— e lPzr Sin ( pz s ) 



(egn 2. 26) 



for Z r > 



And the magnitude of £(^y) i 



is 



| <T)| = 



- Vr~ 



sin(Pz 5 ) 



(eqn 2. 27) 



For the case where Z^ > ly > it can be shown tnat 



q(y) = - 



7 r 
/ 4 



V27T 



VT //0* 



£ 



S//7 






(eqn 2.23) 



and the magnitude of Qi'y) is 



9 ( 7 ) 



rrr* 

= — Sin(pz r ) 



(eqn 2.29) 



A close look at the magnitude cf the pressure spectrum func- 
tion in equation 2.27 reveals that nodes occur for values of 



6Z s = nTT , n- 0 , 1 , 2 , 3 ,... 



(eqn 2.30) 
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as a 



If the magnitude of the pressure spectrum is plotted 
function of the vertical wavenumber then the spacing between 
nodes is uniform and the depth of the source can be 

derived from the relationship given by 



_ 7r 

Zs ~ A P (egn 2. 31) 

In the event the source is deeper than the receiver, looking 
at equation 2.29, we can see that nodes now occur for values 
of 



8Z r n tt , n = 0,1,2, 3,... 



(eqn 2.32) 



Plotting equation 2.29 as a function 
umber, the spacing between nodes 
receiver depth based on 



of the vertical vaven- 
wi-11 now reveal the 



Z 



r 



TT 

A6 



{ egn 2. 33) 



Notice in equations 2.31 and 2.33 the exchange of source and 
receiver depths. By placing the receiver at a shallower 
depth than the source, no new information (namely the source 
depth) is to be found. 

In summary, there is a way to calculate the complex 
pressure amplitude as a function of range using the complex 
pressure spectrum. Conversely, if the complex pressure i?(r) 
has been measured at the receiver, then the complex pressure 
spectrum £ (k) can be derived by making use of the relation- 
ship between the Fourier-Bess el Transform Pairs (see equa- 
tion 2.12), and the depth of the source can be found from 
the magnitude of that spectrum (see equation 2.3 1). 
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The foregoing is the mathematics involved in deriving 
the theoretical pressure spectrum. One purpose of this 
thesis was to compare the theoretical value of the pressure 
spectrum as derived with the Fourier-Bessel Transform with 
the pressure spectrum obtained by use of the computerized 
Discrete Fourier Transform (DFT) , otherwise known as the 
Fast Fourier Transform, or FFT. At this point, it will be 
helpful to review what happens in the FFT. 

Starting with the Fourier- Bessel transform pair shown in 
equation 2.12, f (r) and ( §( / y) can be re-written in terms 
of a discrete sum as is done in the DFT: 
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(egn 2. 34) 



where N ~ 1,2,3,..., up to some large positive integer, and 
represents the number of transform points, and 



yr = (n Ay) (m A r) = 



2jTirm 

N 



(egn 2. 35) 



based on the relationship between A r and Ay 
in sampling theory [ Bef . 2: p. 2], 

Ay Ar = 

or 

. 27T 

= ~NA? 



as described 



(egn 2.36) 



By convention, the sign of the exponential function in equa- 
tion 2.34 is taken as negative when performing the "forward" 
transform. 
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(egn 2.37) 
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and is taken as positive when performing the "inverse" 
trans form. 



FFT 



-I 



g(n AyJ 



= £(m&r) 



(egn 2. 38) 
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D. THE EFFECTS OF SUBF1CE ROUGHNESS 

Surface roughness [Ref. 5: p.409] reduces the Lloyds 

Mirror effect. The rougher the surface, the greater the 
effect for angles increasingly closer to grazing incidence. 
If H represents the surface roughness factor, then equation 

2.3 can be re-written as 

ikR 



P(R) = 



m 2 ikR 
-M e 



R 



R 



where 
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(egn 2. 40) 
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Cp is the grazing angle of incidence, H is the average 

height of the surface gravity waves, and \ is the 

wavelength of the narrowband continuous wave acoustic 
signal. When M is less than 1, the surface is said to be 
smooth; when M is greater than 1, the surface is said to be 
rough . 

looking at Figure 2.1, we ran see that the overall 
effect of increasing surface roughness is to reduce the 
contribution of the surface reflection or "image” to the 
interference pattern at the receiver. In other words, as H 
increases, the energy detected by the receiver is increas- 
ingly only the energy coming directly from the source. When 
the effect of surface roughness is a maximum such that the 
contribution to the pressure field of the "image" is 
completely suppressed, then 

f°° JkR, 

£s (y) = Jo Jo(ir)rcjr 



e iQ\*s-*r\ 

~ 



(egn 2.41) 



and the magnitude of the spectrum is 



I sCy) i = | 



(egn 2. 42) 



E. THE EFFECT OF ADDING NOISE 

In an effort to simulate a realistic ocean environment, 
a normally distributed narrowband noise field based cn a 
Rician distribution [Bef. 6: p. 189] was adapted. The 
resulting acoustic pressure field can be expressed as 
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^(R) = £(R) + n 



(egn 2. 43) 



where n, is a normally-distributed random function with a 
mean cf zero, a standard deviation of one and possesses both 
amplitude and phase. The noise operates independently on 
both the amplitude and phase components of the pressure 
field. Since n is independent of range and wavenumber, it 
is treated as a constant by the Fourier Transform. 

A fourth intention of this thesis was to discover the 
degree to which ncise degrades the wavenumber spectrum 
compared with the corresponding degradation of the range 
domain. To produce the random noise fields used in this 
research, the routine listed as "GGNML" in the IMSL library 
[Ref. 9] was called upon twice to generate independent, 
pseudo-random functions which interact separately and simul- 
taneously with both the amplitude and phase components of 
?{R) . Uhile not an elaborate scheme, owing to time and 

resource constraints, it was felt that this modest simu- 
lacrum of narrowband noise would give a fairly accurate, 
first cut ’’feel" for the effects of noise on the wavenumber 
spectrum. The reader should be aware, however, that no 
direct consideration cf coherency along the range path was 
taken into account by this method. 
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Figure 2.1 The Geometry of the Lloyds Mirror Effect 
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Figure 2.3 The Relationship Between K, y and Q 
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III. RESULTS AND CONCLUSIONS 
A. TBE FFT ALGORITHB 

The research for this thesis was conducted entirely with 
computer algorithms to model the Lloyds Mirror, the Fast 
Fourier Transform, the range windows, waves, and the narrow- 
band noise. This was necessary since the Lloyds Mirror is 
an idealized representation of a situation that occurs only 
rarely in nature, and at that is limited to isospeed surface 
ducts. Also, time and money constraints were such that the 
use of computer models was an absolute necessity. For these 
reasons and because the FFT is the heart of the Wavenumber 
Technique, finding a reliable, easy-to-use computer algo- 
rithm to perform the FFT on the complex pressure function 
was deemed very important. 

Initially, the IMSL library routines, FFT2C and FFTCC 
[Ref. 9: p. 232], were used to generate the pressure 
spectra. Although the author was unable to ascertain just 
how Stamey performed the FFT on his data, it is highly prob- 
able that he used an IMSL library routine. While in nis 
thesis, Blanchard stated that he used the FFTCC routine. 
The magnitude of a typical spectrum generated by the routine 
FFT2C is shown in Figure 3. 1. 

The same scenario used to generate Figure 3. 1 was used 
to produce the spectrum shown in Figure 3.2, but the routine 
FFTCC was utilized. As can be seen, the resulting graphs of 
the pressure spectrum magnitude are practically identical. 
Indeed, the set-up of the pressure fields for insertion into 
each routine differ only in that FFT2C requires a data set 
consisting of an integral power of two number of points, 
whereas FFTCC will handle any number of points. And, when 
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the Cooley-Tukey FFT listed in Appendix I was used to trans- 
form the same pressure field as was utilized to produce 
Figures 3.1 and 3.2, Figure 3.3 was the result. This algo- 
rithm differs somewhat from the IMSL routines by the manner 
in which the pressure field data is treated. Both FFT2C and 
FFTCC are designed to handle a complex array. The 
Cooley-Tukey FFT, on the other hand, is designed to trans- 
form real data; but the pressure field is complex. To get 
around this apparent conundrum, the pressure data is defined 
as either real or imaginary by its placement in the array. 
As can be seen, the resultant graph of the spectrum magni- 
tude is almost indistinguishable from that produced by 
either IMSL routine. In conclusion, it does not matter 
which FFT routine is used to generate the pressure spectrum, 
as long as the pressure field in the range domain is 
"fitted" into the array in a manner suitable to the 
particular algorithm used. 



TABLE I 

Critical Values Osed in the Research 



K (m -1 ) 


N 


r (m) 


(m 3 ) 


3.0 


1024 


0.324 


1.172 x 10~ 2 


2.0 


1024 


0.785 


7.816 x 10" 3 


1.0 


1024 


1.571 


3.906 x 10 -3 



It is important to note four things about the FFT at 
this point. 

1. Since the FFT is given only a finite sample of a 
function while the theoretical Fourier Transform 
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looks at the complete function, the computer algo- 
rithm is designed to overcome this necessary short- 
coming by assuming the submitted finite sample 
repeats its pattern an infinite number of times. 
This can present problems. If the right hand side of 
one pattern does not flow smoothly into the left hand 
side of the next repetition, as is shown in Figure 
3.4, then high frequency oscillations known as the 
Gibbs Phenomenon can be introduced into the spectrum. 
(See Figure 3.5) This is a very real danger where 
acoustic pressure fields are concerned. For example: 
if the sample size does not include the entire 
pattern (i.e. , the pressure magnitude has not been 
completely attenuated prior to the endpoints of the 
sample) , this unwanted "jitter" in the spectrum can 
result. To avoid this situation, not only was a 
pressure field symmetric about tne origin (the 
source) used, but also a Hanning Window was 
constructed to reduce both the amplitude and the 
phase components of the pressure field to zero at the 
endpoints, thus avoiding the Gibbs Phenomenon. 
Figure 3. 6 shows the same pressure field depicted in 
Figure 3.4 combined with a Hanning Window. Its spec- 
trum is shown in Figure 3.7. The Gibbs Phenomenon is 
absent. The form of the window is given below. The 
reader’s attention is directed to Appendix I for a 
view of the manner in which it is combined with the 
pressure field. 

The Hanning Window: 

W • 0.s\ I + COS<Mg.) (egn 3.1) 
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where ’’NPTS" is the number of data points, 
and I = 1 ,2,3, .. . ,NPTS. 

The Hanning Window does not affect the acoustic pres- 
sure field in any other way; the interference pattern 
in the range domain characteristic of the Lloyds 
Mirror for an Acoustic Doublet remains the same 
except for its amplitude near the end points. As 
shall be shown later in this chapter, it is the 

pattern of the entire pressure field, not simply the 
amplitude of the pattern, on which the R.T. performs 
its legerdemain. 

2. Because the vertical wavenumber beta is real and has 

physical meaning for values of the horizontal 

wavenumber less than k, only those values of the 
magnitude of the acoustic pressure spectrum —corre- 
sponding to < k were retained. 

3. In the Wavenumber Technigue, the relationship between 

the horizontal range step size /\r and the hori- 
zontal wavenumber step size /\ r y is extremely 

important. Secalling equation 2.36, 



(egn 3.2) 

for a given "N, 11 the term on the right hand side, 
2 7T* /N , is a constant. Consequently, 




Ar « (£xy) 



(egn 3.3) 



Simultaneously, in accordance with sampling theory, 
in order to prevent aliasing in the pressure spec- 
trum, there is an upper limit on the size of /\r: 
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(egn 3.4) 



A r 



X 



2 



Normally, LA r is chosen such that it is much less 
than half a wavelength; specifically, < X/5 or 

X/6, thus a good sample of the range domain is 
ensured. However, from equation 2.36, the smaller 
^r becomes, the larger /\^y becomes and the 
coarser the sample sizes of the spectrum become. 
Also, to ensure adequate samples in the spectrum, 
examination of equations 2.25 and 2.36 and [Bef. 10] 
reveals that 



at - M¥)‘ 



(eqn 3.5) 



where Z is the deeper of the two depths, whether 
source or receiver. So there is a delicate balance 
to maintain between the sizes of /Xr and Ay. in 
the spirit of compromise, the investigator chose to 
use a range step size equal to a quarter of a wavel- 
ength where 



A 



(eqn 3.6) 



For most cases investigated, values for k , N , /\r 

and At are summarized in Table I. As can be seen, 
using these values and equation 3.5, for a " Z 11 value 
of 25 wavelengths, the maximum size At can assume 
and still ensure an adequate sample size is 
At = 0.0004 /m. Clearly the sizes of At shown 
in Table I are too coarse by a factor of approxi- 
mately 20. Only by increasing N , the number of 
sample points, can this discrepancy be corrected. 
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However, as mentioned earlier, computer resources 
were limited- The effect on the research results is 
not significant; nevertheless, it is recommended the 
reader bear this limitation in mind during the 
succeeding pages. 

4. Looking again at Figures 3.1 through 3.3, the FFT is 
the curve marked by the solid line, while the theo- 
retical Fourier Transform is the curve marked by the 
dotted line. Notice the large discrepancy between 
the theoretical values and the actual values in the 
region corresponding to small vertical wavenumbers. 
This happens to correspond to a region of the spec- 
trum characterized by more rapid fluctuations in the 
spectrum. However, the size of beta is constant. 
Consequently, fewer samples of this region of the 
spectrum are available as compared with that region 
corresponding to high beta values. Since the plot- 
ting routine utilized is interpolating .between 
sampled points of the spectrum, the graph at this end' 
is less smooth and results in a marked difference 
between the curves. Where the curves are determined 
by more points and are smoother, the discrepancy 
between theory and computerized reality is much less 
noticeable. 

B. SOURCE DEPTH DETEE3 INATION 



Theory predicts that, for a source shallower than the 
receiver, when the magnitude of the pressure spectrum is 
plotted as a function of the vertical wavenumber, the resul- 
tant nodal spacing can be used to determine source depth by 
either consulting a chart (see Figure 3.8) or by performing 
a simple calculation (see equation 2.31). To test this 
prediction, several scenarios were designed. Figures 3.9 
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through 3.11 are typical of the many cases run. In each. 



source depth is less than 


receiver depth. 


For each case. 


careful measurement of the 


A/3 spacing coupled with use of 


equation 2.31 


yielded the 


(air eady 


known) 


source depth with 


less than a two per cent 


error. 


which can be 


attributed to 


human error in 


the measurement of 


A P- 










TABLE II 








Eesul 


ts of Source 


Depth Determination 


2 uns 


Figure Nr. 


ZS (m) 


ZR (M) 


(Me 


as ’d) 


Z (Calc 1 d) 


3.9 


21.998 


47.124 


0.14 




22.4 m 


3.10 


31.416 


47.124 


0.10 


1 f 


31.4 " 


3 . 11 


15.708 


47.124 


0.20 


1 f 


15.7 " 


3.12 


47.124 


21.998 


0.14 


1 1 


22.4 " 


3 . 13 


47 . 124 


31.416 


0.10 


I 1 


31.4 " 


3 . 14 


47 . 124 


15.708 


0.20 


1 1 


15.7 " 



To test what happens when the receiver depth is shal- 
lower than the source depth, the cases depicted in Figures 
3.12 through 3.14 were run. The scenario is the same as for 
Figures 3.9 through 3.11; only the source and receiver 
depths have been exchanged. As can be seen, the A Q 
spacing corresponds with the receiver depth. This is again 
as predicted by theory. Table II summarizes the findings of 
Figures 3.9 through 3.14. 
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The scenarios shown in Figures 3.9 through 3.14 all 
start with a pressure field measured horizontally from the 
source; i.e., the exact range to the source is known. What 
if the range to the source is not known, or only minimal 
information concerning the range is known? In this case, 
only a portion of the pressure field starting at seme 
initial range E e can be sampled. What is the effect on 
the spectrum? 

To answer these questions, a range window was 
constructed and placed at varying initial ranges from the 
source. As Figure 2.2 illustrates, as range from the source 
increases, the nodes of the Lloyds Mirror interference 
pattern in the range domain become more widely spaced. 
Given fewer nodes to sample, will the spectrum be the same? 
In accordance with theory, the spectrum is independent of 
range (see equation 2.39) and so should be the same wherever 
the window is placed. 

Figures 3.15 through 3.17 illustrate that the spectrum 
is not entirely independent of range. The figures utilize 
the same scenario as was used above, differing only in the 
range from the source at which the sampling of the pressure 
field begins. As can be seen, beginning at the right hand 
side which corresponds to high vertical wavenumber values, 
or the low horizontal wavenumber values (see equation 2.8), 
a "washing out" of the spectrum occurs, increasing as the 
range window is moved further from the source. Again, the 
theoretical curve is marked by the dotted line. 

When the magnitude of the pressure spectrum is plotted 
as a function of the vertical wavenumber, the useful portion 
of the resulting graph is the last two-thirds of the beta 
range. For K = 2.0/m, this is the range 



fully in Section III A 4. Kith this caveat in mind and 



0.67/m < 




This phenomenon is explained more 



based upon multiple source, receiver, range window 
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combinations run through the model, the following criterion 
was established: the nodal spacing of the spectral plots 

was no longer deter cinable once the initial range H 0 was 
increased to approximately three times the source depth. In 
other words, minimum knowledge of the range from receiver to 
source must be available to ensure adequate sampling of the 
acoustic pressure field in the range domain. This, in turn, 
will produce a pressure spectrum of such a quality that 
source depth determination can be made. 

What causes this not entirely unexpected result? As the 
range window is moved further from the source, for 

Zj ^ Z P , those samples of the pressure field corresponding 
to low gamma values are lost to the spectrum first. Low 

gamma values correspond to waves arriving at high angles 
with respect to the horizontal. At greater receiver ranges, 
5 more closely approximates r and the arrival angles become 
closer tc the horizontal (see Figure 2.1). 

As the window is moved further from the source and fewer 
nodes of the interference pattern are encountered, if a 
wider window were used, could this "washing out" effect in 
the spectrum be minimized or even eliminated? In the 
author’s opinion, it would be minimized; because of the 
reasoning in the preceding paragraph, it probably would not 
be eliminated. However, a wider window would require more 
sample points; because of equation 2.36, one cannot merely 



to this investigator were limited, and the test cases run 
were already using the maximum available array sizes, thus 
precluding a quantitative illustration of this hypothesis. 
Succeeding research into the W.T. might include ways to test 
this point. 

In conclusion, the Wavenumber Technique has the poten- 
tial for being a valuable operational tool in determining 
the depth of an acoustic source provided the receiver is at 



make 




Time and the computer resources available 



a greater depth than the source. It is strongly recommended 
that further research into use of the W.T. in a realistic 
ocean environment be done. 



C. THE EFFECT OF SURFACE ROUGHNESS 

This and the succeeding section were purely qualitative 
portions of the research. Therefore, in the figures the 
spectra were plotted as functions of the horizontal 
wavenumber for ease of discussion except where a particular 
point about the vertical wavenumber was being illustrated, 
and the theoretical curve was omitted. 

As described in Chapter II, the effect of surface rough- 
ness is to suppress the contribution of the image to the 
dipole interference pattern. Moreover, by increasing H, the 
effect is to suppress the image’s contribution to the pres- 
sure field in the range domain and reduce the pressure spec- 
trum to the contribution of the direct path wave only (see 
equation 2.41); the resultant magnitude of the spectrum is 
inversely proportional to the vertical wavenumber (see equa- 
tion 2.42). As Figures 3.13 through 3.21 illustrate, this 
is exactly what happens in the FFT. In Figure 3.21, an 
inverse beta curve has been manually superimposed onto the 
actual curve to reflect this point. 

Looking at the pressure field in the range domain jP(R) 
the effect of surface roughness is inversely proportional to 
the wavelength of the signal. For longer wavelengths, the 
Lloyds Mirror effect is more tolerant of surface roughness 
than it is for shorter wavelengths. Does the same relation- 
ship hold for the pressure spectrum? Comparison of the 
respective sea states in Figures 3.22 and 3.23 with Figures 
3.24 and 3.25 illustrate that it does. 

This investigator also looked at the effect of varying 
the range window in the presence of waves. This is similar 
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to what was looked at in Section B of this chapter. Results 
should be similar and for the same reasons. For a given sea 
state, source and receiver geometry and wavelength, the 
range window was moved successively further from the source 
(see Figures 3.26 through 3.29). As can be seen, the 
’’washing out” effect of the spectral nodes begins to effect 
the lower horizontal wavenumber values first. 

D. TBE EFFECT OF BOISE 

Random noise was simulated in this research. Being 
independent of range and frequency, it was treated as a 
constant by both the Fourier Transform and the FFT. The 
investigator wished to compare qualitatively the effects of 
equal amounts of noise on the pressure field in the range 
domain and in the wavenumber domain. To do this, succes- 
sively larger values of the "noise factor" were used in 
computing n (see equation 2.43). Values of the noise 

''“v/ 

factor used in the research were limited to the maximum 
amplitude cf the pressure field £(£) • Since dipole radia- 
tion is predicated, as range from the source increases, the 
maxima of the interference pattern decrease as the inverse 
square of the range (i.e., R ^). Hence, even small noise 
factors can have a devastating effect on the pressure field 
as the receiver is stepped out in range. The destructive 
effect on the spectrum of successively more intense amounts 
of noise at varying wavelengths is illustrated in Figures 
3.30 through 3.47 

Several conclusions concerning the pressure spectra can 
be drawn from these results: 

1. For equal amounts of noise, longer wavelengths are 
affected less than shorter wavelengths. Compare 
Figures 3.31 and 3.32 with Figures 3.43 and 3.44, and 
Figures 3.34 and 3.35 with Figures 3.46 and 3.47. 
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This was cot a surprising result since the sane prin- 
ciple holds for the range domain. 

The spectrum is affected less by noise than is the 



Figures 3.36 through 3.47 for an illustration of this 
point at two different wavelengths and various noise 
levels. This was a rather welcome surprise. 

3. If the magnitude of the pressure spectrum is plotted 



be detected even after the pressure function in the 
range domain has been "swallowed up" by the noise 
field. Compare Figures 3.33 through 3.35 and Figures 
3.42 through 3.44. This was not surprising m view 
of 2. above. 



E* SUMMATION 

All of the foregoing results are based specifically on 
the lloyds Mirror for an Acoustic Doublet. This is a highly 
idealized and artificial environment. Any thought of imme- 
diately applying conclusions reached in this paper to a real 
world situation is premature. However, the results are- 
still significant if only for supporting the statement made 
by each preceding investigator that the Wavenumber Technique 
should be examined very closely for a future operational 
role, especially in view of the current trends in source 
levels and ambient ncise levels. 

Certain fundamental conclusions regarding the author's 
research into the Wavenumber Technique can now be made 
within the confines of the above restriction: 

1. Source depth can be determined quickly and easily 
from the acoustic pressure spectrum provided 
a) the receiver is deeper than the source. 



range domain. Compare Figures 3.30 through 3.35 and 



as a function of beta, the 




spacing can still 



2 . 



b) some knowledge of the range from receiver to 
source is available, and 

c) the magnitude of the pressure spectrum is plotted 
as a function of the vertical wavenumber. 

Source depth determination in even a noisy environ- 
ment is possible. While the introduction of 
increasing amounts of narrow band noise adversely 
affects both the range domain and the wavenumber 
domain, the pressure spectrum seems able to withstand 
the chaotic destruction longer than does the pressure 
field in the range domain. 

3. The ability to determine source depth is adversely 
affected by the height of surface gravity waves since 
surface roughness reduces the Lloyds Mirror as the 
sea state increases. 
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PRESSURE SPECTRUM VS KZ 

SOURCE DEPTH = 20 m, K = 2.0/m, 
RECEIVER DEPTH = 50 m, N = 1024 




Figure 3. 1 Pressure Spectrum Osing FFT2C 
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PRESSURE SPECTRUM MAGNITUDE 
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PRESSURE SPECTRUM VS KZ 

SOURCE DEPTH = 20 m, K = 2.0/m, 
RECEIVER DEPTH = 50 m, N = 1024 




Figure 3.2 Pressure Spectrum Using FFTCC 
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PRESSURE SPECTRUM MRGNITUDE 
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PRESSURE SPECTRUM VS KZ 



SOURCE DEPTH = 20 m, K = 2.0/m, 
RECEIVER DEPTH = 50 m, N = 1024 



o 




Figure 3.3 Pressure Spectrum Using Cooley-Tukey FFT 
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Figure 3.4 Hon-smoothed Pressure Field 
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PRESSURE SPECTRUM MAGNITUDE 
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Figure 3.5 Pressure Spectrum Showing Gibbs Phenomenon 
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Figure 3.6 Pressure Field Combined With a Hanning Window 
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PRESSURE SPECTRUM VS KH 
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HORIZONTRL WRVENUMBER, 1/M 



Figure 3.7 Pressure Spectrum Combined With a Hanning Window 
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Figure 3.8 Theoretical Source Depth Determination Curve 
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PRESSURE SPECTRUM MAGNITUDE 
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PRESSURE SPECTRUM VS KZ 



RECEIVER DEPTH = 47.124 M, K = 2.0/M, 
N = 1024 




VERTICAL WAVENUMBER, 1/M 



Fiqure 3.9 Graph of Pressure Spectrum, 
Source at 22.0 Meters 
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PRESSURE SPECTRUM MAGNITUDE 



PRESSURE SPECTRUM VS KZ 

RECEIVER DEPTH = 47.124 M, K = 2.0/M, 
N = 1024 




Figure 3.10 Graph of Pressure Spectrum, 
Source at 31.4 Meters 
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PRESSURE SPECTRUM VS KZ 

RECEIVER DEPTH = 47.124 M, K = 2.0/M, 
N = 1024 
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Figure 3.11 Graph of Pressure Spectrum, 
Source at 15.7 Meters 
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PRESSURE SPECTRUM MAGNITUDE 



PRESSURE SPECTRUM VS KZ 

SOURCE DEPTH = 47.124 M, K = 2.0/M, 
N = 1024 
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Figure 3-12 Graph of Pressure Spectrum, 
Receiver at 22.0 Meters 
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PRESSURE SPECTRUM VS KZ 

SOURCE DEPTH = 47.124 M, K = 2.0/M, 
N = 1024 
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VERTICRL WRVENUMBER, 1/M 



Figure 3.13 Graph of Pressure Spectrum, 
Receiver at 31.4 Meters 
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PRESSURE SPECTRUM VS KZ 

SOURCE DEPTH = 47.124 M, K = 2.0/M, 
N = 1024 
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VERTICRL WRVENUMBER, 1/M 



Figure 3.14 Graph of Pressure Spectrua, 
Receiver at 15.7 Meters 
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PRESSURE SPECTRUM VS KZ 

SOURCE DEPTH = 21.998 M, K = 2.0/M 
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Figure 3. 15 Pressure Spectrum, 
Range Hindow Set at 47.1 Meters 
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PRESSURE SPECTRUM VS KZ 

SOURCE DEPTH = 21.998 M, K = 2.0/M, 
Receiver Depth = 47.124 M, N = 1024 



o 




Figure 3- 16 Pressure Spectrum, 
Range Window Set at 50.3 Meters 
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PRESSURE SPECTRUM MUGNITU0E 



PRESSURE SPECTRUM VS KZ 

SOURCE DEPTH = 21.998 M, K = 2.0/M, 
RECEIVER DEPTH = 47.124 M, N = 1024 



o 




Figure 3. 17 Pressure Spectrum, 
Range Window Set at 62.8 Meters 
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PRESSURE SPECTRUM MAGNITUDE 



1 



PRESSURE SPECTRUM VS KZ 

SOURCE DEPTH = 21.998 M, K = 2.0/M 
RECEIVER DEPTH = 47.124 M, H = 0.0 M 




VERTICAL WAVENUMBER, 1/M 



Figure 3.18 Pressure Spectrum, Sea State 0 
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PRESSURE SPECTRUM MAGNITUDE 



PRESSURE SPECTRUM VS KZ 



SOURCE DEPTH = 21.998 M, K = 2.0/M 
RECEIVER DEPTH = 47.124 M, H = 1.0 M 




VERTICRL WAVENUMBER, 1/M 



Figure 3.19 Pressure Spectrum, Sea State 2 
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PRESSURE SPECTRUM VS KZ 



SOURCE DEPTH = 21.998 M, K = 2.0/M 
RECEIVER DEPTH = 47.124 M, H = 3 M 



o 




VERTICRL WAVENUMBER, 1/M 



Figure 3.20 Pressure Spectrum, Sea State 3 



65 



PRESSURE SPECTRUM MAGNITUDE 



PRESSURE SPECTRUM VS KZ 



SOURCE DEPTH = 21.998 M, K = 2.0/M 
RECEIVER DEPTH = 47.124 M, H = 10 M 




VERTICRL NRVENUMBER, 1/M 
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Figure 3.21 Pressure Spectrum, Sea State 5 
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PRESSURE SPECTRUM MAGNITUDE 
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PRESSURE SPECTRUM VS KH 

SOURCE DEPTH = 21.998 M, N = 1024 
j RECEIVER DEPTH = 47.124 M 
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Figure 3,22 Pressure Spectrum, Sea State 0, K = 1.0 
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PRESSURE SPECTRUM MRGNITUDE 



i 

Figure 



PRESSURE SPECTRUM VS KH 

SOURCE DEPTH = 21.998 M, N = 1024 
RECEIVER DEPTH = 47.124 M 




3.23 Pressure Spectrum, Sea State 3, K = -1.0 
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PRESSURE SPECTRUM MRGN1TUDE 



PRESSURE SPECTRUM VS KH 

SOURCE DEPTH = 21.998 M, N = 1024 
RECEIVER DEPTH = 47.124 M 




Figure 3.24 



Pressure Spectrum, Sea State 0, K = 2.0 



PRESSURE SPECTRUM VS KH 

SOURCE DEPTH = 21.998 M, N = 1024 
RECEIVER DEPTH = 47.124 M 




Figure 3.25 Pressure Spectrum, Sea State 3, K = 2.0 
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PRESSURE SPECTRUM MAGNITUDE 



1 



. PRESSURE SPECTRUM VS KH 

SOURCE DEPTH = 21.998 M, N = 1024 
RECEIVER DEPTH = 47.124 M 




Figure 3.26 Pressure Spectrum, 
SS 2, Range Window Set at 50 Meters 
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PRESSURE SPECTRUM MRGNITUDE 



PRESSURE SPECTRUM VS KH 



SOURCE DEPTH = 21.998 M, N = 1024 
RECEIVER DEPTH = 47.124 M 



o 




Figure 3.27 Pressure Spectrum, 

SS 2 , E'ange Bindow Set at 100 Meters 
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PRESSURE SPECTRUM MRGNITUDE 



i 



PRESSURE SPECTRUM VS KH 

SOURCE DEPTH = 21.998 M, N = 1024 
RECEIVER DEPTH = 47.124 M 




Figure 3.28 Pressure Spectrum, 

SS 2, Range Window Set at 200 Meters 
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PRESSURE SPECTRUM VS KH 

SOURCE DEPTH = 21.998 M, N = 1024 
RECEIVER DEPTH = 47.124 M 
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Figure 3.29 Pressure Spectrum, 

SS 2, Range Hindow Set at 300 Meters 
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MAGNITUDE OF PRESSURE AS A FN OF RANGE 

SOURCE DEPTH = 21.998 M, RECEIVER DEPTH = 47.124 M, 
RANGE STEP SIZE = 1.57 M, N = 1024 

CD 
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Figure 3.30 Pressure Field, K = 1.0, jJL = 0.005 
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PRESSURE SPECTRUM MRGNITUDE 
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PRESSURE SPECTRUM VS KH 

SOURCE DEPTH = 21.998 M, K = l.O/M 
RECEIVER DEPTH = 47.124 M, u = 0.005 



o 




H0RIZ0NTRL WAVENUMBER, 1/M 
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Figure 3.31 Pressure Spectrum vs. Gamma 
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PRESSURE SPECTRUM VS KZ 

SOURCE DEPTH - 21.998 M, K = 1 O/M 
RECEIVER DEPTH = 47.124 M, y = 0.005 




VERTICAL WAVENUMBER, 1/M 



Figure 3-32 Pressure Spectrum vs. Beta 
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MRGNITUDE OF PRESSURE RS R FN OF RRNGE 

SOURCE DEPTH = 21.998 M, RECEIVER DEPTH = 47.124 M, 
RANGE STEP SIZE = 1.57 M, N = 1024 
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1—= j 



Figure 3.33 Pressure Field, K = 1.0, fJL = 0.01 
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PRESSURE SPECTRUM MAGNITUDE 



PRESSURE SPECTRUM VS KH 



SOURCE DEPTH = 21.998 M, K = l.O/M 
RECEIVER DEPTH = 47.124 M , u = O.Ol 



o 




Figure 3.34 Pressure Spectrum vs. Ga«ma 
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PRESSURE SPECTRUM MAGNITUDE 






PRESSURE SPECTRUM VS KZ 

SOURCE DEPTH = 21.998 M, K = l.O/M 
RECEIVER DEPTH = 47.124 M, v = O.Ol 



o 




Figure 3.35 Pressure Spectrum vs. Beta 
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MAGNITUDE OP PRESSURE, P(R) 
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MAGNITUDE OF PRESSURE AS A FN OF RANGE 

SOURCE DEPTH = 21.998 M, RECEIVER DEPTH = 47.124 M 
RANGE STEP SIZE = 0.785 M, N = 1024 
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Figure 3.36 



Pressure Field, K = 2.0, 0.0001 



PRESSURE SPECTRUM MRGNITUDE 



PRESSURE SPECTRUM VS KH 

SOURCE DEPTH = 21.998 M, K = 2.0/M 
RECEIVER DEPTH = 47.124 M, u = O.OOOl 




Figure 3-37 Pressure Spectrum vs. Gamma 



82 



PRESSURE SPECTRUM MAGNITUDE 



PRESSURE SPECTRUM- VS KZ 

SOURCE DEPTH = 21.998 M, K = 2.0/M 
RECEIVER DEPTH = 47.124 M, u = O.OOOl 




VERTICAL WAVENUMBER, 1/M 



Figure 3.38 Pressure Spectrum vs. Beta 



83 



r 



i 



MAGNITUDE OF PRESSURE RS R FN OF RRNGE 

SOURCE DEPTH = 21.998 M, RECEIVER DEPTH = 47.124 M, 
RANGE STEP SIZE = 0.785 M, N = 1024 

LD 




Figure 3.39 Pressure Field, K = 2.0, £4= 0.001 
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PRESSURE SPECTRUM MRGNITUOE 




Figure 3.40 Pressure Spectrum vs. Gamma 
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PRESSURE SPECTRUM VS KZ 

SOURCE DEPTH = 21.998 M, K = 2.0/M 
RECEIVER DEPTH = 47.124 M, u = 0.001 
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Figure 3.41 pressure Spectrum vs. Beta 
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MAGNITUDE OF PRESSURE RS R FN OF RANGE 

SOURCE DEPTH = 21.998 M, RECEIVER DEPTH = 47 124 M 
RANGE STEP SIZE =0.785 M, N = 1024 
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Figure 3.42 Pressure Field, K = 2.0, fj, = 0.005 
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PRESSURE SPECTRUM VS KH 

SOURCE DEPTH = 21.998 M, K = 2.0/M 
RECEIVER DEPTH = 47.124 M, y = 0.005 
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Figure 3.43 Pressure Spectrum vs. Gamma 
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PRESSURE SPECTRUM MAGNITUDE 



PRESSURE SPECTRUM VS KZ 

SOURCE DEPTH = 21.998 M, K = 2.0/M 
RECEIVER DEPTH = 47.124 M, u =0.005 
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VERTICAL WAVENUMBER, 1/M 



Pigure 3.44 Pressure Spectrum vs. Beta 
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MAGNITUDE OE PRESSURE, P(R] 



MAGNITUDE OP PRESSURE AS A EN OE RANGE 

SOURCE DEPTH = 21.998 H, RECEIVER DEPTH = 47.124 M, 
RANGE STEP SIZE = 0.785 M, N = 1024 
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Figure 3.45 Pressure Field, K = 2. 0, fX~ 0.01 
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PRESSURE SPECTRUM MAGNITUDE 



PRESSURE SPECTRUM VS KH 

SOURCE DEPTH = 21.998 M, K = 2 O/M 
RECEIVER DEPTH = 47.124 M, u = 0.01 
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Figure 3-46 Pressure Spectrum vs. Gamma 
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PRESSURE SPECTRUM MRGNITUDE 



r 



PRESSURE SPECTRUM VS KZ 

SOURCE DEPTH = 21.998 M, K = 2.0/M 
RECEIVER DEPTH = 47.124 M , y = O.Ol 




Figure 3-47 Pressure Spectrum vs. Beta 
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